Electron spin relaxation in rippled graphene with low mobilities 
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We investigate spin relaxation in rippled graphene where curvature induces a Zeeman-like spin- 
orbit coupling with opposite effective magnetic fields along the graphene plane in K and K' valleys. 
The joint effect of this Zeeman-like spin-orbit coupling and the intervalley electron-optical phonon 
scattering opens a spin relaxation channel, which manifests itself in low-mobility samples with the 
electron mean free path being smaller than the ripple size. Due to this spin relaxation channel, with 
the increase of temperature, the relaxation time for spins perpendicular to the effective magnetic 
field first decreases and then increases, with a minimum of several hundred picoseconds around 
room temperature. However, the spin relaxation along the effective magnetic field is determined by 
the curvature-induced Rashba-type spin-orbit coupling, leading to a temperature-insensitive spin 
relaxation time of the order of microseconds. Therefore, the in-plane spin relaxation in low-mobility 
rippled graphene is anisotropic. Nevertheless, in the presence of a small perpendicular magnetic 
field, as usually applied in the Hanle spin precession measurement, the anisotropy of spin relaxation 
is strongly suppressed. 

PACS numbers: 81.05.ue, 85.75.-d, 71.70.Ej, 75.70.Tj 



I. INTRODUCTION 



In recent years, graphene has attracted much inter- 
est due to its potential for the all-carbon based elec- 
tronics and spintronics. 1-14 A number of experiments 
on spin relaxation in graphene on SiC>2 substrate have 
been carried out, with spin relaxation times t s of the 
order of 10-100 ps reportedi^iiSr— Some works sug- 
gested the Elliott- Yafet (EY)^ 3 - mechanism to be dom- 
inant in spin relaxation , 15 i 17 ~ — while the surface chemi- 
cal doping experiment supported the importance of the 
D'yakonov-Perel' (DP) 24 one. 22 These experiments have 
triggered intensive theoretical studies on spin relaxation 
in graphene.— ~— With inversion asymmetry, possibly 
caused by a perpendicular electric field or curvature, 
a Rashba-type spin-orbit coupling (SOC), 34 which cou- 
ples spin to pseudospin, arises in graphenej 29 ' 35 " — It was 
found that with this SOC, the DP mechanism dominates 
the electron spin relaxation. 27 Besides, it was revealed 
that for the EY mechanism, t s /t p ~ n e , where t p is 
the momentum relaxation time and n e is the electron 
densityj 27 i 30 This makes the assessment, which attributes 
the observed linear relation between spin relaxation time 
r s and momentum relaxation time r p with the increase 
of electron density n e to the EY mechanism, in the ex- 
perimental work o 15 i 17 ~ — inappropriate M- The magnitude 
of the Rashba-type SOC caused by a moderate perpen- 
dicular electric fiel d 38 i 39 or curvatur o 29 i 35 is of the or- 
der of fj,eV and the corresponding t s limited by the DP 
mechanism is as large as microseconds^ The adatoms 
were suggested to enhance the local SOC to the order 
of 10 meV and hence provide a possible origin of the 
observed short spin relaxation time . 25 ' 26 i 32 ' 33 i 39 ' 40 Re- 
cently, we set up a random Rashba model incorporat- 
ing the effect of adatoms 3 - 3 - and fitted the experimental 



data from different groupsj 18 ' 19 : 21 : 22 We suggested that 
the DP mechanism dominates spin relaxation in graphene 
and can result in either linear or inversely linear relation 
between r, and t„ 
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Very recently, Jeong et al. reported that curvature 
in graphene can lead to not only the Rashba-type SOC 
which is off-diagonal in the pseudospin space, but also an 
additional SOC diagonal in the pseudospin spaced 9 - This 
additional SOC serves as a Zeeman-like term with oppo- 
site magnetic fields in K and K' valleys, similar to the 
case in carbon nanotube.— ~— Starting from the effective 
SOC, Jeong et al. studied spin relaxation in a chemically- 
clean corrugated graphene where the electron mean free 
path is much larger than the spatial range of the random 
spin-orbit fluctuation and the spatially averaged SOC is 
zero. Under such condition, the spin relaxation is solely 
limited by the spin-flip scattering due to the fluctuation 
of the SOC) 31 ' 44 i 45 as in the rippled graphene studied 
by Dugaev et al. where only the Rashba-type SOC was 
considered. 31 It is noted that in the two valleys, even 
though the effective magnetic fields from the Zeeman- 
like term are opposite, their contributions to the spin 
relaxation via the fluctuation-induced spin-flip scatter- 
ing are independent and identical. Besides, it was found 
that the fluctuations of the Rashba-type SOC and the 
Zeeman-like term make comparable contribution to the 
spin relaxation, leading to a spin relaxation time at least 
of the order of 10 nsi^ 9 . 

Nevertheless, for low-mobility samples where the elec- 
tron mean free path is smaller than the spatial range 
of the random spin-orbit fluctuation, the spin relaxation 
is determined by the local curvature-induced SOC4^ For 
each isolated valley, the Rashba-type SOC, together with 
the Zeeman-like term, still leads to a spin relaxation time 
of the order of microseconds in the DP mechanism. How- 



ever, the two valleys are not independent any more due 
to the intervalley scattering. The opposite effective mag- 
netic fields in two valleys, together with the intervalley 
scattering, give rise to another spin relaxation channel. 
This can be understood by a simple model. We label the 
spin vector in each valley as S M , where fi — ±1 stands 
for the valley located at K or K' point. The spin vec- 
tor in each valley precesses around the effective magnetic 
field from the Zeeman-like term with a frequency w^, 
where uj^ = —<*>-^ and = uj. The intervalley scat- 
tering is characterized by a scattering time t v between 
two valleys. The spin vectors S M , with the initial values 
S/x(0) = S-^O) = So, satisfy the rate equations 

S M (t) + S M (t) x uv + [S^t) - S_ M (i)]/r„ = 0. (1) 

When So lies in the plane perpendicular to the effective 
magnetic field, one has the solution 

^S M (i) = 2So e - t/rs (2) 

with r s = 2/(w 2 t„) in the strong intervalley scattering 
limit ojt v « 1. So far this mechanism has not been re- 
vealed in the literature. 

In this work, we study the spin relaxation in the 
low-mobility rippled graphene (the mobility is around 
2 x 10 3 cm 2 /V-s)2ii and take into account the above 
spin relaxation channel. The electron mean free path 
I is smaller than the ripple size £. In the low tempera- 
ture regime where the intervalley electron-phonon scat- 
tering is negligible, the spin relaxation is dominated by 
the Rashba-type SOC and t s is as large as microseconds. 
However, with the increase of temperature, due to the 
above spin relaxation channel, the relaxation time for 
spins polarized perpendicular to the effective magnetic 
field first decreases and then increases, with a minimum 
of several hundred picoseconds around room tempera- 
ture. 

This paper is organized as follows. In Sec. II, we 
present the model and Hamiltonian. In Sec. Ill, we study 
the spin relaxation in rippled graphene based on the ki- 
netic spin Bloch equations (KSBEs)4£~— The effect of 
temperature, impurity density and electron density on 
spin relaxation is investigated. The anisotropy of spin 
relaxation, without and with a small perpendicular mag- 
netic field, is also addressed. We summarize in Sec. IV. 

II. MODEL AND HAMILTONIAN 

The quasi-periodically rippled graphene can be syn- 
thesized by chemical vapor deposition on copper first 
and then transferred to the SiC>2 substrate perpendicu- 
lar to the z-axis3 The surface morphology of the rippled 
graphene is illustrated in the upper panel of Fig. Q] The 
graphene is curved along the x-axis and we define the 
angle between the y-axis and one carbon-carbon bond 
in counterclockwise direction as a (0 < a < 2ir/3). 29 
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FIG. 1: (Color online) Upper panel: schematic of the rip- 
pled graphene curved along the :r-axis. The angle between 
the y-axis and one carbon-carbon bond in counterclockwise 
direction is a (0 < a < 2n/3). Lower panel: schematic of the 
effective magnetic field from the SOC along one circle around 
K and K' points, (a) The effective magnetic field induced by 
the Rashba-type SOC only (£' = 0); (b) The effective mag- 
netic field from both types of the SOC with a — 0. 

The size and height of the ripples are about 50 and I nm 
respectively^ The radius R of curvature is 100-200 nm. 20 
According to Ref. [29|, the local effective SOC induced by 
the curvature reads 

-ffsoc = (k(lit x ®a y ~T y ® <J X ) + fJ,('nI <E> b • er, (3) 

with b = (cos(7r/2 + 3a), sin(7r/2 + 3a), 0) shown in the 
upper panel of Fig. Q] Here t and I are the Pauli matrices 
and unit matrix in the pseudospin space. <r are the Pauli 
matrices in the spin space. The curvature K = R . The 
parameters £ = 0.15 meV-nm and £' = 0.21 meV-nmi 2 ^ 
On the right-hand side of Eq. ([3]), the first term is the 
Rashba-type SOC reported previously 3 - 5 - and the second 
term is the additional Zeeman-like term diagonal in the 
pseudospin spaced 

With the basis laid out in Refs. [25] anddH, the effective 
Hamiltonian of the electron band reads 

H = ^2 [£k<5 SS '+(^ M k+2.9MBlB)- cr s ; 5']c Al ks t C Al k S '+i^mt- 

(4) 

Here c M ks (c M ks^) is the annihilation (creation) operator 
of electrons in \x valley with momentum k (relative to 
the valley center) and spin s (s — ±4). £k = hvfk with 
Vf = 10 6 m/s. g fa 2 is the effective Lande g-factor ) 50 i 51 
fiB is the Bohr magneton and B is an external mag- 
netic field perpendicular to the graphene plane (when it 
is applied, its magnitude is very small and the effect on 
orbital motio n 50 ' 51 is negligible) . The effective magnetic 



3 



field from the SOC is given by 

S^k = C/c(-sin0 k ,cos0 k ,O) + fiC'nh, (5) 

where 6* k is the polar angle of momentum k. The 
second term on the right-hand side of the equation 
plays the role of the Zeeman-like term with the effec- 
tive static magnetic field along b (— b) in K (K') val- 
ley. The effective magnetic fields, without and with 
the second term, are schematically plotted in the lower 
panel of Fig. [TJ where a is set as zero and hence 
b is along the y-aixs. The Hamiltonian H- mt con- 
sists of both the intravalley and intervalley scatterings. 
The former include the electron-impurity,— electron- 
remote-interfacial phonon,— electron-acoustic phonon, 54 
electron-r-E2 ff optical phonon^ as well as electron- 
electron Coulomb scatterings.— The latter include the 
electron-K-A^ optical phonon^ and electron-electron 
Coulomb scatterings^ 



III. SPIN RELAXATION IN RIPPLED 
GRAPHENE 

The KSBEs 4 ^— are utilized to study the spin relax- 
ation in grahene ; 28 ' 33 



dtp^{t) = 9 ty o M k(i)|coh + dt/Vk(*)|s 



(0) 



where p Mk (i) represent the density matrices of elec- 
trons with relative momentum k in valley y, at time 
t. The coherent terms read 3 t p M k(*)|coh = — l[(^k + 
ig/isB) • er,/j M k(i)], where the Hartree-Fock term from 
the Coulomb interaction is neglected due to the small 
spin polarizationi 28 i 33 i 47 The concrete expressions of the 
scattering terms St/Opk^lscat are given in Ref. [28l . By 
solving the KSBEs, one obtains the spin relaxation time 
r s along direction n from the time evolution of spin po- 
larization P(t) — X}uk r 'M/ /ik(*) £r ' n l- br our calcu- 
lation, unless otherwise specified, the initial spin polar- 
ization is P(0) = 10 %, the spin-polarization direction 
n is along the z-axis, the curvature n = 0.01 nm _1 , the 
electron density n e = 7 x 10 11 cm -2 , the impurity den- 
sity rii = 2 x 10 12 cm -2 and the external magnetic field 
B = 0. For spin relaxation along the z-axis, the direction 
of b, determined by angle a, is irrelevant. 



A. Temperature dependence of spin relaxation 

We first study the temperature dependence of the spin 
relaxation. In Fig. [2] the spin relaxation time r s is 
plotted against temperature T at different curvatures 
n. The electron mean free path I is around 25 nm (the 
corresponding mobility is around 2.6 x 10 3 cm 2 /V-s) in 
the whole temperature regime investigated, as shown in 
Fig. [2] with the scale on the right-hand side of the frame. 
Therefore the electron mean free path is always smaller 
than the ripple size. It is indicated by this figure that 



when T < 100 K, r s is of the order of microseconds. How- 
ever, when T goes beyond 100 K, r s decreases rapidly to 
several hundred picoseconds at T ~ 200-300 K. Never- 
theless, when T further increases and exceeds the room 
temperature, r s begins to increase with T. It is also 
seen from the figure that r s reaches its minimum at a 
higher temperature with larger curvature k. Moreover, 
when T = 50 K and T > 250 K, r s is proportional to 
k~ 2 . However, in the intermediate temperature regime 
100 < T < 250 K, the spin relaxation times are nearly 
the same for different values of k. This scenario is un- 
derstood as follows. 
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FIG. 2: (Color online) Solid, dashed and dotted curves: tem- 
perature dependence of spin relaxation time with different 
values of n. Chain curve with the scale on the right-hand side 
of the frame: temperature dependence of mean free path. 



= 7 x 10 11 cm"' and n, = 2x 10 1 



We first focus on the case with k — 0.01 nm -1 . When 
T < 100 K, the intervalley scattering is negligible 2 ^ and 
spins in two valleys relax independently. The spin re- 
laxation is then determined by the Rashba-type SOC. 
The Zccman-likc term only serves as an in-plane effective 
magnetic field which mixes the in-plane and out-of-plane 
spin relaxations^ It is known that with the Rashba- 
type SOC only and strong electron-impurity scattering, 
tx = T||/2 = h 2 /(4( 2 K 2 Tp), where t±_ [t\\) is the out-of- 
plane (in-plane) spin relaxation timej 25 ' 32 Therefore in 
the presence of the effective magnetic field, the spin re- 



laxation time along the z-axis is t s 
h 2 /(3C 2 K 2 T„). With r„ 



2/(r7 



'-L ^ 'II > ~ 

1 /(3( 2 k 2 t p ). With t p = l/vf, one can estimate t s 
2.4 /is at T = 50 K, as shown in the figure. When T 
increases, the intervalley electron-K-A^ optical phonon 
scattering becomes important and opens another spin 
relaxation channel together with the opposite effective 
magnetic fields in two valleys. According to the model 
presented in the introduction, in the weak intervalley 
scattering limit ujt v > 1 (uj = 2Q'n/h for the concrete 
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situation here), one has the solution 



: sin 



(Vw 2 -T„- 2 i-r» (7) 



10' 



with </) = arctan ^w 2 t 2 — 1 when So is perpendicular 
to the effective magnetic field. This indicates that the 
spin relaxation time is solely determined by the inter- 
valley scattering time, i.e., t s — t v . Therefore, in the 
weak intervalley scattering limit, r s decreases with the 
enhancement of the intervalley scattering and hence the 
increase of T. However, in the strong intervalley scat- 
tering limit, t s = 2/(uj 2 T v ) as given by Eq. ([2]) in the 
introduction. In such case, t s increases with the in- 
crease of T. The crossover from the weak to strong in- 
tervalley scattering limit with the increase of T is deter- 
mined by r^ 1 w lu = 2C,'k/H. At this crossover point, 
t s t v w H/(2£'k), which is estimated to be 170 ps, just 
as the value shown in the figure. 

Based on the above analysis, it is understood that in 
the zero intervalley scattering limit r s = h 2 / (AC, 2 k 2 t p ) 
and in the strong intervalley scattering limit r s = 
h 2 1 (2£ /2 k 2 t v ). In both limits r s is proportional to kT 2 . 
However, in the weak intervalley scattering limit with 
tJ 1 -c t^ 1 < u>, t s is determined by the intervalley 
scattering time t v and remains insensitive to k. Besides, 
the crossover point in the nonmonotonic temperature de- 
pendence of spin relaxation time moves to a higher tem- 
perature with the increase of curvature k as determined 
by the relation r" 1 w 2C,'n/h. These properties mani- 
fest themselves in the curves with different values of k in 

Fig.m 



B. Effect of intervalley scattering and SOC on spin 
relaxation 



The intervalley scatterings include both the electron- 
electron Coulomb and electron-phonon scatterings. 2? 
However, the essential role played in the spin relax- 
ation channel revealed in this work is the intervalley 
electron-phonon scattering. That is because the interval- 
ley Coulomb scattering which transfers electrons between 
the two valleys is negligible due to the large momentum 
transfer between them and hence the small scattering 
matrix element. Only the intervalley Coulomb scatter- 
ing which does not lead to any electron transfer between 
the valleys is considered^ 

For comparison, we show the temperature dependence 
of spin relaxation time with different intervalley scatter- 
ings included in Fig. [3] It is seen that the intervalley 
Coulomb scattering is unimportant in the whole temper- 
ature regime under study while the intervalley electron- 
phonon scattering affects spin relaxation effectively. Par- 
ticularly, when the intervalley electron-phonon scattering 
is excluded, t s becomes insensitive to T. That is because 
in the nearly isolated valleys, r s = h 2 /(3C, 2 K 2 T p ) with t v 
being dominated by the electron-impurity scatteringi 28 i 33 
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FIG. 3: (Color online) Temperature dependence of spin relax- 
ation time with the inclusion of different intervalley scatter- 
ings. Solid curve: with both the intervalley electron-electron 
(e-e) and electron-phonon (e-p) scatterings; Dashed curve: 
without the intervalley e-e scattering; Chain curve: without 
the intervalley e-p scattering; Dotted curve: without both 
the intervalley e-e and e-p scatterings. n e — 7 x 10 11 cm -2 , 
m = 2 x 10 12 cm" 2 and « = 0.01 nm" 1 . 
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FIG. 4: (Color online) Temperature dependence of spin relax- 
ation time for the genuine case ^ 0), the case without the 
Zeeman-like term (£' = 0) and the case without the Rashba- 



type SOC (( = 0). n e 
and k — 0.01 nm" . 
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In Fig. [4l we further compare the genuine spin re- 
laxation (CC' 7^ 0) to the ones without the Zeeman- 
like term ((' = 0) and without the Rashba-type SOC 
(C = 0), in order to reveal the effect of the two types 
of SOC on spin relaxation. When the Zeeman-like term 
is absent ((' = 0), the two valleys are degenerate and 
t s = t ± = h 2 /{A( 2 K 2 T p )^^ The dotted curve in Fig. H 
satisfies this relation very well. When the Rashba-type 
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SOC is absent (£ = 0), the spin relaxation is caused 
by the opposite effective magnetic fields jointly with the 
scattering between two valleys. At low temperature, t s 
approaches infinity due to the suppression of intervalley 
scattering. By comparing the two curves with £ = and 
Cf = to the genuine one with ((' ^ 0, one finds that 
in the genuine situation the spin relaxation is dominated 
by the Rashba-type SOC when T < 100 K and by the 
Zeeman-like term when T > 100 K. 



C. Spin relaxation with different impurity and 
electron densities 

We calculate the spin relaxation with a higher impurity 
density rii = 3 x 10 12 cm -2 to explore the effect of impu- 
rity density on spin relaxation. With this impurity den- 
sity the electron mean free path / is decreased (compared 
to the values shown in Fig. [2]) and the condition Z ^ is 
satisfied. In Fig. the ratios of the momentum scatter- 
ing rate and spin relaxation time with rii = 3 x 10 12 cm" 2 
to those with rii — 2 x 10 12 cm" 2 , labeled as ^(t' 1 ) and 
t](t s ), are plotted against the temperature. It is shown 
that ^{t^ 1 ) remains around 1.5 in the whole temperature 
regime under study, as the electron-impurity scattering 
dominates t p . t](t s ) is about 1.5 at T = 50 K and rapidly 
decreases to 1 with the increase of T. That is because 
the spin relaxation is sensitive to the intravalley scatter- 
ing only when t" 1 < tJ 1 (T < 150 K). Particularly, at 
T = 50 K, r s cx r" 1 and rj(r s ) = ^(t" 1 ) = 1.5. When 
T > 150 K, the spin relaxation becomes insensitive to the 
intravalley scattering and hence the increase of impurity 
density. 
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FIG. 5: (Color online) Temperature dependence of the ratios 
of the momentum scattering rate (dashed curve) and the spin 
relaxation time (solid curve) with m = 3 x 10 12 cm -2 to 
those with m — 2 x 10 12 cm -2 , labeled as ^(t^ 1 ) and r/(T B ) 
respectively. n e — 7 x 10 11 cm -2 and n = 0.01 nm -1 . 
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FIG. 6: (Color online) Temperature dependence of spin re- 
laxation time with different electron densities. Solid curve: 
n e = 7 x 10 11 cm" 2 ; Dashed curve: n e = 1.4 x 10 12 cm" 2 . 
m = 3 x 10 12 cm" 2 and k = 0.01 nm" 1 . 



While changing the impurity density influences the 
spin relaxation only at low temperature (T < 150 K), 
the change of electron density is expected to affect spin 
relaxation in a large temperature regime, as both the in- 
travalley and intervalley scatterings depend on the elec- 
tron density. In fact, with the increase of electron den- 
sity, the intravalley electron-impurity scattering is weak- 
ened while the intervalley electron-phonon scattering is 
strengthened. Therefore, with the increase of electron 
density, as long as the condition I -C £ is satisfied, the 
spin relaxation time t s decreases in the weak intervalley 
scattering limit (r s cx t~ 1 when T < 100 K and t s = t v 
when 100 < T < 250 K) but increases in the strong in- 
tervalley scattering limit (r s cx r" 1 when T > 250 K). 
We calculate the spin relaxation with a higher electron 
density n e = 1.4 x 10 12 cm" 2 and the impurity density 
rii = 3 x 10 12 cm" 2 . For such case the electron mean 
free path I has the largest value of 30 nm at T = 50 K. 
In Fig. [SI we compare the spin relaxation for this case to 
the one with n e = 7 x 10 11 cm" 2 and = 3 x 10 12 cm" 2 . 
It is indeed shown that with the increase of electron den- 
sity, t s decreases when T < 250 K but increases when 
T >250 K. 

When the electron density is further increased and 
eventually I > £, the mechanism discussed in this work 
is not valid any more and the spin-flip scattering in- 
duced by the fluctuation of the SOC contributes to spin 
relaxation i 31 i 33 ' 44 i 45 When the latter is dominant, the 
spin relaxation time t s is of the order of 10 ns and insen- 
sitive to the temperature i 29 ' 31 
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D. Anisotropy of spin relaxation without and with 
small perpendicular magnetic field 

The above studies are limited to the spin relaxation 
along the z-axis. In fact, due to the in-plane effective 
magnetic field along b, the spin relaxation time along 
any direction perpendicular to b is identical. Also, ac- 
cording to the simple model presented in the introduc- 
tion, there is no spin relaxation along b. However, due 
to the Rashba-type SOC, spins relax along b exponen- 
tially, with the spin relaxation time t s = ti (ti can be 
obtained from t± given by the dotted curve in Fig. 31 
as T|| = 2tj_). This spin relaxation time is of the order 
of microseconds and insensitive to temperature. Due to 
the difference in spin relaxations along and perpendicu- 
lar to b, the spin relaxation in the plane parallel to b, 
e.g., the graphene plane, is anisotropic. When the initial 
spin-polarization direction n deviates from b by an an- 
gle #n,b, the spin polarization along n relaxes in the form 



P(t) = P(0)[cos 2 



,be 



-t/r. 



: sin 



scattering limit 



_L 



< r. 



l+sm^ n)b /(i)]. Here/(t) = 
in the weak intervalley 

- 1 <uj (100 < T < 250 K) and 



Ty 2 t + 



/(<) = e u r "*/ 2 in the strong intervalley scattering limit 
t" 1 > U) (T > 250 K), according to Eqs. Q and © 
respectively. 



10' 



10" 



10 4 



10' 



10" 











B=10 mT: 6 n b =0 








— e— 




B=0: e n b =0 


- • - 




i 


— X- 


n e =7xlO U cm" 2 
n i= 2xl0 12 cm" 2 
K=0.01 nm" 1 







50 100 150 200 250 300 350 400 
T(K) 



Usually, the spin relaxation in graphene is ex- 
perimentally studied by the Hanle spin precession 
measurement J^r— In such measurement, a small mag- 
netic field B (at most of the order of 100 mT) perpendic- 
ular to the graphene plane is applied. Under the perpen- 
dicular magnetic field, the spin relaxations along different 
directions in the graphene plane are mixed. In the rippled 
graphene studied here, the spin relaxations along both b 
and the direction perpendicular to b are exponential in 
the strong intervalley scattering limit (T > 250 K), with 



the two spin relaxation rates being 1 



and 



r„/2, 



spectively. Therefore, when B > H\t^ 1 — io 2 T v /2\/{2g^B) 
(about 10 mT around room temperature), the spin relax- 
ation along any direction in the graphene plane has the 



unique rate t~ 



= 



II 



uj 2 t v /2)/2 sa UJ 2 T V /A. 5b How- 



ever, this is not the case in the weak intervalley scattering 
limit (T < 250 K), as there is spin oscillation along the 
direction perpendicular to b [refer to Eq. ([7}]. We apply 
a magnetic field of magnitude of 10 mT perpendicular 
to the graphene plane and calculate the spin relaxations 
along b (# n .b = 0) and the direction perpendicular to 
b in the graphene plane (# n .b = 7r /2) respectively at 
different temperatures. The temperature dependence of 
the spin relaxation time along these two directions, with 
and without the perpendicular magnetic field, is plotted 
in Fig. [7J For comparison, we also plot the tempera- 
ture dependence of the parameter, r° , calculated from 
the average value of the perpendicular-magnetic-field- 
free spin relaxation rates along these two directions. It 
is shown that due to the perpendicular magnetic field, 
the spin relaxation becomes isotropic in the graphene 



plane when T > 250 K, with t s 



-t v /A. However, 



when T < 250 K, the anisotropy is strongly suppressed 
but not completely eliminated. Consequently, despite 
the angle a and the measured spin-polarization direc- 
tion, the observed spin relaxation time in low-mobility 
rippled graphene by means of the Hanle spin precession 
measurement is expected to have a marked nonmonotonic 
dependence on temperature, with a minimum of several 
hundred picoseconds located around room temperature. 



IV. CONCLUSION AND DISCUSSION 



FIG. 7: (Color online) Temperature dependence of spin re- 
laxation time along different directions in the graphene plane 
with (solid curves) and without (dashed curves) an external 
magnetic field along the z-axis. Solid curve with squares 
(open circles): spin relaxation along the direction parallel 
(perpendicular) to b with external magnetic field of 10 mT 
along the z-axis; Dashed curve with closed circles (triangles): 
spin relaxation along the direction parallel (perpendicular) to 
b without external magnetic field. The dotted curve with 
crosses stands for the temperature dependence of the pa- 
rameter, t°, calculated from the average value of the spin 
relaxation rates along directions parallel and perpendicular 
to b without external magnetic field. n e = 7 x 10 11 cm -2 , 
n,=2x 10 12 cm" 2 and /t = 0.01 nm" 1 . 



In conclusion, we have studied the electron spin relax- 
ation in rippled graphene with low mobilities. The elec- 
tron mean free path is smaller than the ripple size and 
the spin relaxation is determined by the local SOC in- 
duced by curvature. The curvature not only leads to the 
Rashba-type SOC, but also supplies a Zeeman-like term 
with opposite effective magnetic fields along graphene 
plane in two valleys^ We show that this Zeeman-like 
term, together with the intervalley electron-K-A^ optical 
phonon scattering, gives rise to a spin relaxation channel 
in rippled graphene at high temperature (with tempera- 
ture T > 100 K). This spin relaxation channel can cause 
a marked nonmonotonic dependence of spin relaxation 
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on temperature along the direction perpendicular to the 
effective magnetic field from the Zeeman-like term. A 
minimal spin relaxation time of several hundred picosec- 
onds is obtained around room temperature. However, the 
Rashba-type SOC dominates the spin relaxation along 
the effective magnetic field, leading to a temperature- 
insensitive spin relaxation time of the order of microsec- 
onds. Therefore, anisotropy exists in the spin relaxation 
in the graphene plane. In the presence of a small per- 
pendicular magnetic field in Hanle spin precession ex- 
periment, the anisotropy of in-plane spin relaxation is 
strongly suppressed. Particularly, in the strong interval- 
ley scattering limit, the spin relaxation in the graphene 
plane becomes isotropic, with the spin relaxation time 
being two times as large as that along the direction per- 
pendicular to the effective magnetic field in the absence 
of the perpendicular magnetic field. 

The spin relaxation channel revealed in this work man- 
ifests itself only when the electron mean free path is 
smaller than the ripple size. However, it is noted that 
when this spin relaxation channel is dominant (at tem- 
perature T > 100 K), the spin relaxation is insensitive 
to the intravalley scattering, mainly contributed by the 
electron-impurity scattering in the low-mobility sample. 
It is also noted that in reality the ripples in graphene may 
not be quasi-periodic as proposed in this work. The radii 
of the curvatures may be spatial dependent^ Therefore, 
the observed spin relaxation is coherently averaged over 
the curvature morphology in the low-mobility sampled 
In such case, the spin relaxation is determined by the 
mean square of curvature in the zero and strong inter- 
valley scattering limits. The spin relaxation channel re- 
vealed in this work should also exist in carbon nanotudc. 



However, in the work of Semenov et al. where the spin 
relaxation in carbon nanotube was studied^ the inter- 
valley scattering was not considered and this spin relax- 
ation channel was not incorporated. 

Finally, we discuss the experimentally observed spin 
relaxation in the low-mobility rippled graphene by Avsar 
et al.££ In their sample with an electron density of 
7.5 x 10 11 cm -2 , the in-plane spin relaxation time ob- 
tained from the Hanle spin precession measurement in- 
creases mildly from about 130 to 150 ps when T increases 
from 5 to 300 Ki 2 ^ Avsar et al. estimated the effect of 
the local SOC from curvature and excluded it from the 
spin relaxation mechanism (refer to the supplementary 
information of Ref . [20h . However, it is noticed that they 
only took into account the Rashba-type SOC according 
to Ref. HH and the resulting spin relaxation time is natu- 
rally much longer than their experimental data. The spin 
relaxation channel revealed in this work supplies a pos- 
sible origin of the experimentally observed short spin re- 
laxation time around room temperature. To account for 
the experimental data in the whole temperature regime, 
some other factors, e.g., the randomly enhanced SOC by 
the substrate and/or adatoniS ) 25 ' 26 ' 32 ' 33 ' 39 ' 40 may need to 
be considered. 
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